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Abstract 

We present the complete set of Renormalisation Group Equations (RGEs) at one loop for the 
non-exotic minimal U{1) extension of the Standard Model (SM). It includes all models that are 
anomaly-free with the SM fermion content augmented by one Right-Handed (RH) neutrino per 
generation. We then pursue the numerical study of the pure B — L model, deriving the triviality 
and vacuum stability bounds on an enlarged scalar sector comprising one additional Higgs singlet 
field with respect to the SM. 
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I. INTRODUCTION 



The Large Hadron Collider has essentially been built to confirm or disprove the existence 
of one or more Higgs bosons. A lot of effort has therefore been put into studying models 
that can accommodate the Higgs mechanism of Electroweak Symmetry Breaking (EWSB), of 
which the Higgs (pseudo) scalar particles are remnants. Among these, the most studied one 
is the Standard Model (SM). Unfortunately, the SM is flawed. There is now experimental 
evidence of new phenomena that cannot be explained by the SM, notably (very small) 
neutrino masses. At the same time, it should be noted that the accidental U{1)b-l global 
symmetry (where B{L) is the Barion(Lepton) number) is not anomalous in the SM with 
massless neutrinos, but its origin is not understood. It thus becomes appealing to extend the 
SM to explain simultaneously the existence of both neutrino masses and the B — L global 
symmetry by gauging the U{1)b-l group and subject it to spontaneous EWSB induced 
by the Higgs mechanism (therefore generating a massive Z' state on the same footing as 
massive and Z states are generated from the breaking of the SU{2)l x f/(l)y gauge 
symmetry of the SM). Consequently, this requires that the fermion and scalar spectra are 
enlarged to account for gauge anomaly cancellations in such a way as to evade direct searches. 
Minimally, this requires the addition of a scalar singlet and three (massive) right-handed 
neutrinos, one per generation l|-l4|, the latter entering the see-saw mechanism to explain 
the smallness of the detected (SM-like) neutrino masses 



a. 



Generally, U{1) gauge factors mix and the mixing is controlled by further gauge couplings 
6|. These extra parameters can then be reabsorbed in an effective parameterisation [TJ. In 



" the SM, spanning over several 
" fl fl 

^, |8[ is a particularly simple 



the model we are considering, two Abelian groups are present, and just one extra coupling is 
therefore needed to account for the mixing, which can effectively be reabsorbed. The arising 
model is a minimal (i.e., one-dimensional) U{1) extension o 
benchmark models among which the "pure" B — L model 
example, as we will describe in the next Section. 

Unsurprisingly, because of the simple nature of such a model, following the experimental 
results on neutrino masses, a plethora of papers have been published studying the phe- 
nomenology of the B — L model at colliders. They have dealt with the detectabilty of the 
Z' boson (see jlj] for earlier studies on generic hadron colliders) at the LHC 8-12| and at 
a future Linear Collider (LC) 13l-ll5|. some analyses concentrating on the Z' decaying via 
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heavy neutrinos, in particular into three [8| and four [16| leptons in the final state, with 



distinctive displaced vertices due to long lived neutrinos, a clear signature of physics beyond 
the SM. Also, the testability at the LHC of the see-saw mechanism in this model has been 
evaluated in detail fl?! . 

In comparison to the gauge and neutrino sectors, the Higgs part of this model has un- 
dergone much less scrutiny. Apart from the benchmark study of 9(] dating back a few years, 
only recently a systematic analysis of the Higgs sector of the B — L model has started, in 



the attempt to define the boundaries of the associated parameter space. Ref. [18| dealt with 
the limits stemming from the imposition of perturbative unitarity on the model. Here, we 
intend to pursue further into this attempt, by investigating the triviality and vacuum stabil- 
ity conditions ensuing in the B — L model, through a RGE analysis aiming at defining the 
physical values of the masses and couplings of the two Higgs states emerging in the model 
after EWSB, the latter depending upon the maximum energy scale after which also such a 
scenario ceases to be valid and further new physics dynamics ought to be invoked. We believe 
that, with the LHC now on line, it is of paramount importance to theoretically constrain 
the Higgs sector of a new physics scenario that, while incorporating the SM, it remedies its 
major fiaw without leading to a proliferation of new particles and/or interactions, thereby 
retaining much of the predictivity and testability of the SM. 

In the past and yet recent years, a lot of effort has been spent for similar studies. For 



reviews on the SM and on some of its extensions, see Refs. 



19 



20| and references therein. 
Concerning the study presented here, earlier works focusing on extra singlet scalars or Eq- 



insp ired U{1) augmented gauge groups in non-suspersimmetric 21 



23 



22l | and suspersimmetric 



24| extensions of the SM, respectively, have already been considered. 
The plan of the paper is as follows. In section [TTl we describe the model under study. In 
section UTTl we describe our computational techniques. The following section |IV] presents our 
numerical results while we conclude in section |Vl We also have an appendix, where we list 
the RGEs of the model that we have dealt with. 

II. THE PARAMETERISATION 

n 

The model under study is the minimal U{1)b-l extension of the SM (see ref. [4] for 
conventions and references), in which the SM gauge group is augmented by a f/(l) factor. 
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related to the Baryon minus Lepton {B — L) gauged number. In the complete model, the 
classical gauge invariant Lagrangian, obeying the SU{3)c x SU(2)l x f/(l)y x U{1)b~l 
gauge symmetry, can be decomposed as: 

^ = + .^YM + + -^Y ■ (1) 

The scalar Lagrangian is: 

= {D^H)^ D,H + {D^x)^ D,x - V{H, x) , (2) 
with the scalar potential given by 



V{H,x) = m^H^H + ix'lx? + [WH \ x? 
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'■H^H + f,' I X r +Xi{H^Hy + X, I X r +hH^H \ x T , (3) 



where H and x s^i'^ the complex scalar Higgs doublet and singlet fields, respectively. 

We generalise the SM discussion of spontaneous EWSB to the more complicated classical 
potential of eq. (j3]). To determine the condition for V{H,x) to be bounded from below, it 
is sufficient to study its behaviour for large field values, controlled by the matrix in the first 
line of eq. ([3]). Requiring such a matrix to be positive-definite, we obtain the conditions: 

4AiA2-A^>0, (4) 

Ai,A2>0. (5) 

If the above conditions are satisfied, we can proceed to the minimisation of V as a function 
of constant Vacuum Expectation Values (VEVs) for the two Higgs fields. Making use of 
gauge invariance, it is not restrictive to assume: 

W-Q). (X)-^. (6) 

with V and x real and non-negative. The physically most interesting solutions to the min- 
imisation of eq. ([3]) are obtained for v and x both non- vanishing: 



To compute the scalar masses, we must expand the potential in eq. ([3]) around the minima 
in eqs. ([7]) and ([8]). We denote by hi and h2 the scalar fields of definite masses, nih^ and 
rufi^ respectively, and we conventionally choose m^^ < m1^. After standard manipulations, 
the explicit expressions for the scalar mass eigenvalues and eigenvectors are: 



m 



^ = \y + Aax' - v/(Ait;2 - A2x2)2 + [X^xvf , (9) 



2 - \y + X^x" + ^{\iv^ - A2x2)2 + {X^xvY , (10) 

(11) 





H 







cos a — sm « 
sin a cos a 




where -f < a < f fulfils^: 

sin 2a = ^ ^'^"^ , (12) 

v/(Ait;2 - A2x2)2 + (A3a;t;)2 

„ Aif ^ - Aax^ 

cos/a = — , (Id) 
^{Xiv^ - X2X^y + [X^xvY 

For our numerical study of the extended Higgs sector, it is useful to invert eqs. (jH]), ( 1TU|) 

and ( |T2l) . to extract the parameters in the Lagrangian in terms of the physical quantities 

rrihj^, and sin 2a: 

2 2 

Ai = -4(l-cos2a) + -4(l + cos2a), 

2 2 
TTl TTT 

A2 = ^(l-cos2a) + ^(l + cos2a), 

^---(^)- (") 

Moving to the J^ym, the non-Abelian field strengths therein are the same as in the SM 
whereas the Abelian ones can be written as follows: 



^YM = --F^'F.u - ^F^^'F' , (15) 



where 



Ff^u = d,,B, - d,B^ , (16) 
F' = d,Bl - d^B' . (17) 



^ In all generality, the whole interval < a < 27r is halved because an orthogonal transformation is invariant 
under a a + w. We could re-halve the interval by noting that it is invariant also under a — — a if we 
permit the eigenvalues inversion, but this is forbidden by our convention m^^ < mf^^. Thus a and —a are 
independent solutions. 
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In this field basis, the covariant derivative is: 

D, = d, + igsT'^G^^ + igT-W; + igiYB, + i(gY + g[YB-L)B' . 



(18) 



To determine the gauge boson spectrum, we have to expand the scalar kinetic terms as 
for the SM. We expect that there exists a massless gauge boson, the photon, whilst the 
other gauge bosons become massive. The extension we are studying is in the Abelian sector 
of the SM gauge group, so that the charged gauge bosons will have masses given by 
their SM expressions, being related to the SU{2)l factor only. Using the unitary-gauge 
parameterisation, the kinetic terms in eq. ([2]) become: 



{D''H)^D^H = ]-d^hd^h + \{h + vf{Ql)\gW,''aa + giB^' + gB'^'^ 
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M|2 



and 



+ {gW/-g,B^^-gB'>^f 



(19) 



(20) 



where we have taken K 



B-L 



2 in order to guarantee the gauge invariance of the Yukawa 
terms (see eq. (123]) ). In eq. ffT^ we can recognise immediately the SM charged gauge bosons 
W^, with Mw = gv/2 as in the SM. The other gauge boson masses are not so simple to 
identify, because of mixing. In fact, in analogy with the SM, the fields of definite mass are 
linear combinations of B'^, and B'^. The explicit expressions are: 

— sin cos t?' sin sin t?' ^ ^ 
cos cos ^' — cos 'dw sin ^' 
sin -t}' cos 

















sin i?, 




\b'^) 






with -f < ^9' < f , such that: 



(21) 



tan 2^9' 



2g^/? 



9l 



g2 + lQ[^)2g'2_g2_g2 



(22) 



and 



Ma 



9' + gl 



1 /fi ^ + 16(f)V 
9^ + 9l 



sin 2^9 + 9l 



(23) 
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where 

. 2g,/^^T7i 

sm 2v = — . 

^ {g' + 16(f -g'- glY + {2gY{9^ + qD 
LEP experiments 25j constrain < 10"'^. Present constraints on the VEV x (see sec- 
tion IIVBI) allow a generous range of Tj. 

The fermionic Lagrangian (where k is the generation index) is given by 

3 

fc=i 

MlkLl^D'^lkL + lekRltiD^GkR + iVkRlf.D^'vkR^ , (24) 

where the fields' charges are the usual SM and B — L ones (in particular, B — L = 1/3 
for quarks and —1 for leptons with no distinction between generations, hence ensuring 
universality). The B — L charge assignments of the fields as well as the introduction of 
new fermionic RH heavy neutrinos (i^ij's) and a scalar Higgs field (x, charged +2 under 
B — L) are designed to eliminate the triangular B — L gauge anomalies and to ensure the 
gauge invariance of the theory, respectively. Therefore, the B — L gauge extension of the 
SM gauge group broken at the Electro- Weak (EW) scale does necessarily require at least 
one new scalar field and three new fermionic fields which are charged with respect to the 
B — L group. 

Finally, the Yukawa interactions are: 

■^Y = -yjkqJIdkRH - yJ^qJZUkRH - y]j^ljLekRH 



-VjJjLJ^kRH - Vji, {lyRYjiykRX + h.c. , (25) 

where H = ia^H* and i,j,k take the values 1 to 3, where the last term is the Majorana 
contribution and the others the usual Dirac ones. 

Neutrino mass eigenstates, obtained after applying the see-saw mechanism, will be called 
ui (with / standing for light) and Uh (with h standing for heavy), where the first ones are 
the SM-like ones. 



A. Realistic models 



The generic model that has been previously introduced spans over a continuous set of 
minimal U{1) extensions of the SM, that can be labelled by the properties of the charge 
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assignments to the particle content. Notice that such models are, by construction, all and 
only those that are not anomalous with the SM fermion content augmented by one RH 
neutrino per generation. Therefore, many ii^g-inspired U{1) scenarios, such as U{1)^ or 
U{l)rj, are not included in our generic model. 

Free parameters in our parameterisation are those appearing in the covariant derivative 
of eq. ( !T8|) . We recall the Abehan part only: 



This form of the covariant derivative can be re-written defining an effective coupling and 
an effective charge qe'- 



As any other parameter in the Lagrangian, ]) and g[ are running parameters [g, |7|, 
therefore their values ought to be defined at some scale. A discrete set of popular Z' models 
(see, e.g., Refs. [2^ can be recovered by a suitable definition of both g and g[. 



We will focus our numerical analysis on the scalar sector of the "pure" B — L model, that 
is defined by the condition g{QEw) = 0, i.e., we nullify it at the EW scale. This implies no 
mixing at the tree-level between the B — L Z' and SM Z gauge bosons. Other benchmark 
models of our general parameterisation are for example the Sequential SM (SSM), defined 
by = Y (that in our notation corresponds to the condition g[ = at the EW scale) 
and the f/(l)K model, for which RH fermion charges vanish (that is recovered here by the 
condition g = —2g[ at the EW scale). 

It is important to note that none of the models described so far is orthogonal to the U{1)y 
of the SM, therefore the RGE running of the fundamental parameters, g and g[, will modify 
the relations above. The only orthogonal U{1) extension of the SM is the "S'O(10)-inspired" 
f/(l)^ model, that in our notation reads g = —^g'l- Although the g and g[ couplings run 
with a different behaviour, the EW relation g/g[ = —4/5 is preserved (at one-loop) at any 
scale. 

Nonetheless, as indeed true for the SM, the gauge sector affects marginally the scalar 
sector in its running, so the analysis we are going to show is effectively independent of the 
specific charge assignation. However, we might expect differences regarding the interplay 
between the gauge sector and the neutrino evolution, that impinge on the vacuum stability 
studies of the scalar sector as the top quark does for the SM Higgs sector. We will report 



D^ = --- + ig^YB^ + i{gY + g[YB-L)B' 



gEY^ = gY + g[YB^L. 



(26) 
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separately on the study of the differences in the RGE study of the specific benchmark models 
in our generic parameterisation. 

III. COMPUTATIONAL DETAILS 

The complete set of RGEs for the generic model are derived for the parameters in the 
Lagrangian and are collected in appendix |A1 For their numerical study, we put boundary con- 
ditions at the EW scale on the physical observables: m/j^, mh^, a, v, Mz', g[,g,ml;^'^, that 
we trade for m, fi, Ai, A2, A3, x, 1/^23 using, for the relevant parameters therein, eq. ffH]) . 
Where stated in the text, we impose boundary conditions on some parameters of the La- 
grangian rather than on the physical observables. This is done for consistency of those 
studies. 

For the pure B — L model, object of the numerical analysis in this work, the definition 

5^ = holds, and as a consequence, we also have that the B — L breaking VEV x can be 

Mz' 

easily related to the new Z' boson mass by x = — where we fixed g'^ = 0.1. Regarding 

2^71 

the neutrinos, for simplicity we consider them degenerate and we fix their masses to m];'^'^ = 
TTiu^ = 200 GeV (whenever not specified otherwise), a value that can lead to some interesting 
phenomenology [_8J. The free parameters in our study are then m^^, rrih^, a and x. The 
general philosophy is to fix in turn some of the free parameters and scan over the other ones, 
individuating the allowed regions fulfilling the following set of conditions. 

We first define a parameter to be "perturbative" for values less than unity. This is a 
conservative definition, as we could relax it by an order of magnitude and still get values 
of the parameters for which the perturbative series will converge^. RGE evolution can then 
constrain the parameter space of the scalar sector in two complementary ways. From one 
side, the couplings must be perturbative. This condition reads: 

< Ai,2,3(Q') < 1 Vg'<Q, (27) 

and it is usually referred to as the "triviality" condition. On the other side, the vacuum of 
the theory must be well-defined at any scale, that is, to guarantee the vahdity of eqs. (jl]) 

^ Notice that, in analogy with QED, the parameters upon which the perturbative expansion is performed 
are usually of the form ^/a = g/\/^, rather then being g itself. 
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and ([5]) at any scale Q' < Q: 

< Ai,2,3(Q') and 4Ai(g')A2(Q') - \1{Q') > V Q' < Q . (28) 

Eq. ( |28|) is usually referred to as the "vacuum stability" condition. In contrast to the SM, 
in which it is sufficient the Higgs self-coupling A be positive, in the case of this model the 
vacuum stability condition (and especially the second part of eq. (128!) ) can be violated even 
for positive Ai,2,3- 

One should notice that our conventional choice nih^ < nih^, as noted previously, let us 
consider a and —a as two independent solutions, although the theory is manifestly invariant 
under the symmetry a — )■ —a. These two solutions are complementary, meaning that the 
region excluded by the choice rrih^ < rrih^ at a certain value of the angle a is precisely the 
allowed one for the complementary angle 7r/2 — a. The special case a = 7r/4 is symmetric, 
and corresponds to maximal mixing between the scalars. a = corresponds to a SM scalar 
sector totally decoupled from the extended one, and hi is the usual SM Higgs boson, a = 7i/2 
is the specular case, in which h2 plays the role of the SM Higgs boson. 

Notice also that, again in contrast to the SM in which the gauge couplings have a marginal 
effect, in our case the RH neutrinos play for the extra scalar singlet the role of the top quark 
for the SM Higgs in the vacuum stability condition^. Their RGE are then controlled by the 
Yukawa coupling with a negative contribution coming from g[ (see eq. ( 1A9I) ). Therefore, in 
some regions of the parameter space, the impact of the gauge sector is not marginal and 
can effectively stabilise the otherwise divergent evolution of the Majorana Yukawa couplings 
for the RH neutrinos. We will report on the effect of RH neutrinos in our analysis in 
section IIVDI 

A final remark is in order about eq. flA33p . the evolution of A3, the mixing parameter 
of the scalar potential (see eq. ([3])). This RGE is almost proportional to A3 itself, so a 
vanishing boundary condition is almost stable^. Non-proportional terms arise from the new 
gauge couplings (g and g[), i.e., deviations from the vanishing boundary conditions are of 
the order of the gauge coupling, hence quite small. They are particularly negligible in the 

^ Also notice that we have three RH neutrinos, as we have three colours for the top quark. However, they 
are Majorana particles rather than Dirac ones, so they carry half (independent) degrees of freedom than 
the top quark. 

^ From the last line of eq. (|T4)) . setting A3 = corresponds to a = 0, but not vice versa. 
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pure B — L model, as also g has a vanishing boundary condition, with a weak departure 
from it due to the mixing in the gauge coupling sector j4|. Nonetheless, other benchmark 
models in our general parameterisation could show different behaviours. 



IV. RESULTS 

We present here our results for the pure B—L model, the chosen benchmark of our general 
parameterisation. We will first present a brief analysis of the gauge sector, followed by a 
quick review of the present experimental constraints on the Higgs boson masses. Finally, we 
will fully describe the scalar sector analysis, argument of this paper^. 

A. Gauge sector 

Before starting the analysis of the scalar sector, we can briefly look at the gauge sector, 
where the RGE evolution gives us indications for the validity of the model concerning the 
gauge couplings. In particular, their evolution must stay perturbative up to some particular 
scale. In the B — L model, the conditions that the free parameters in the gauge sector must 
fulfil are: 

g[{Q') < 1 V Q' < Q and g{QEw) = , (29) 

where the second condition in eq. (129|) defines the pure B — L model. 

Varying the scale Q, the maximum scale up to which we want the model to be well-defined, 
we get an upper bound on g'i{QEw) as a function of Q, as shown in figure [H Typical results 
are summarised in table HI 



LogioiQ /GeY) 


3 


5 


7 


10 


15 


19 


3i{Qew) 


0.860 


0.693 


0.593 


0.497 


0.397 


0.342 



TABLE I: Maximum allowed values by eq. ^2M) for g[{QEw) iri the B — L model for 

selected values of the scale Q. 



^ Notice that we study the gauge sector of the model (and, in particular, the Abelian part of it) indepen- 
dently of any other sector as the corresponding RGEs fully decouple. 
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FIG. 1: Maximum allowed values by eq. (d^j for g[{QEw) the B — L model as a 

function of the scale Q. 



B. Experimental limit 

Past and current experiments have set limits on the scalar sector parameters in the SM 



for LEP and Ref. 



for 



as well as in various extensions of it, see for example Ref. 
Tevatron. For the model discussed here, the relevant analysis is summarised in figure |2l 
in which a generic overall factor ^ has been introduced. Such parameter is defined as the 
coupling(s) to the Z boson of the Higgs particle(s) in the considered extension normalised 
to the SM: 

e-fS^, (30) 

hence it parametrises the deviations of the new model with respect to the SM. 

In the minimal U{1) extension of the SM, argument of this paper, two scalar eigenstates 
exist: the one coming from the Higgs singlet, required to break the extra U{1)b~l gauge 
factor (and therefore giving the Z' gauge boson a mass), and the one coming from the Higgs 
doublet, required to break the SM gauge symmetry to give masses to the W and Z bosons. 
With reference to eq. f|TT]) . we called hi the lightest of such eigenstates, that couples to 
the Z boson proportionally to cos a, and with /i2 we referred to the heaviest scalar, that 
couples to the Z boson proportionally to sin a. Hence, the LEP lower bounds on the scalar 
masses of the U{1)b-l extension here considered are read straightforwardly from figure [2] 
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20 40 60 80 100 120 



mj,(GeV/ci, 

FIG. 2: The 95% C.L. upper bound on ^ = gnzz/gflzz ^J- the B - L model, 

C, = cosQ;(sina) for H = hi{h2). 

by considering: 

£ = cos a ioT H = hi , 

(31) 

^ = sin a for H = h2 , 
i.e., the limit for /ii(/i2) are extracted by considering ^ as the cosine(sine) of the mixing 
angle in the scalar sector (see eq. ( iTTl) and the following ones). 

Figure [2] shows the lower bound on the Higgs mass as a function of ^. The SM Higgs is 
recovered by the condition ^ = 1. We see that we can have significant deviation from the 
SM Higgs mass limit, > 114.4 GeV, only for values of the angle a > 71/4, for the lightest 
state hi. For example, for a = it/ 3, the LEP limit on the lightest Higgs state reads as 
nihj^ > 100 GeV. That is, in this model, a light Higgs with mass smaller than the SM limit 
can exist only if it is highly mixed, i.e., the light Higgs is mostly the singlet state. For the 
same value of the angle, the limit for rrih^ is more stringent than the condition m^^ > rrihi, 
in fact for a = 7r/3, nih,^ ^114 GeV must be fulfilled. 

The LEP experiments are also able to provide a lover bound for the B — L breaking VEV 
X. In fact, the LEP bound on the B — L Z' mass |30|, 

^ > 7 TeV (32) 
can be rewritten clS Qb lower bound for the VEV: 

X > 3.5 TeV , (33) 
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since Mz' = '2xg[ in the pure B — L model. 
C. Scalar sector 

Given the simphcity of the scalar sector in the SM, the triviality and vacuum stability 
conditions can be studied independently and they both constrain the Higgs boson masses, 
providing an upper bound and a lower bound, respectively. In more complicated models as 
the one considered here, it might be more convenient to study the overall effect of eqs. fl27p - 
fl25]) . since there are regions of the parameter space in which the constraints are evaded 
simultaneously. This is the strategy we decided to follow. 

Figure [3] shows the allowed region in the parameter space mh^-nih^ for increasing values of 
the mixing angle a, for fixed VEV x = 7.5 TeV and heavy neutrino masses m,y^ = 200 GeV, 
corresponding to Yukawa couplings whose effect on the RGE running can be considered 
negligible. For a = 0, the allowed values for rrihi SM ones and the extended scalar 

sector is completely decoupled. The allowed space is therefore the simple direct product of 
the two, as we can see in figure [3al When there is no mixing, the bounds we get for the new 
heavy scalar are quite loose, allowing a several TeV range for m/ij, depending on the scale 
of validity of the theory. We observe no significant lower bounds (i.e., rrih,^ > 0.5 GeV), as 
the RH Majorana neutrino Yukawa couplings are negligible. 

As we increase the value for the angle, the allowed space deforms towards smaller values 
of rrih^. If for very small scales Q of validity of the theory such masses have already been 
excluded by LEP, for big enough values of Q, at a small angle as a = 0.1, the presence 
of a heavier boson allows the model to survive up to higher scales for smaller hi masses if 
compared to the SM (in which just hi would exist). Correspondingly, the constraints on 
mfi2 become tighter. Moving to bigger values of the angle, the mixing between hi and /12 
grows up to its maximum, at a = it /A, where hi and /i2 both contain an equal amount of 
doublet and singlet scalars. The situation is therefore perfectly symmetric, as one can see 
from figure [3cl Finally, in figure |3dl we see that the bounds on m/jj are getting tighter, 
approaching the SM ones, and those for are relaxing. That is, for values of the angle 
71 /A < a < 7t/2, the situation is qualitatively not changed, but now /i2 is the SM-like Higgs 
boson. Visually, one can get the allowed regions at a given angle 7r/2 — a by simply taking 
the transposed about the nih-^ = rrih^ line of the plot for the given angle a. 
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FIG. 3: Allowed values in the rrih^ vs. rrih.^ space in the B — L model by eqs. [2l\ ) and ^2^, 
for ( T^al) a = 0, [SU) a = 0.1, ( T^cj) a = tt/4 and [Sdl) a = 7t/3. Colours refer to different 
values of Q / GeV: blue (10^), red (lOi^), green (10^^), purple (10^^) and cyan (10^^). The 
shaded black region is forbidden by our convention rrih^ > irih^ , while the shaded red region 
refers to the values of of the scalar masses forbidden by LEP. Here: x = 7.5 TeV, 

m^^ = 200 GeV. 

Per each value of the angle, we can then fix the lighter Higgs mass m^-^ to some benchmark 
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values (allowed by LEP for the SM Higgs) and plot the allowed mass for the heavier Higgs 
as a function of the scale Q. This is done in figure HI where the allowed masses are those 
contained between the same colour lines. Notice that here the VEV x is fixed to a different 
value, X = 3.5 TeV. The effects of changing the VEV x will be described in section IIVEI 

As previously noticed, the allowed range in ruh^ gets smaller as we increase the angle. 
Apart from the case a = where there is no dependency at all from m^^, there is a strong 
effect from rrihj^ on the bounds on rrih^- Not all the allowed regions at a fixed hi mass are 
contained in the region for a smaller m^^. This is true only for rrih^ > 160 GeV. For smaller 
m/i^'s, the distortion in the allowed region constraints tightly m/^j for the survival of the 
model to big scales Q. This is because such distortion is just towards smaller hi masses, see 
figure [31 

Complementary to the previous study, we can now fix the light Higgs mass at specific, 
experimentally interesting^, values, i.e., ruh-^ = 100, 120, 160 and 180 GeV, and show the 
allowed region in the rrih^ vs. a plane. This is done in figure [H 

From this figures it is clear the transition of h2 from the new extra scalar to the SM-like 
Higgs boson as we scan on the angle. As we increase (up to = 160 GeV), a bigger 
region in rrih^ is allowed for the model to be valid up to the Plank scale (the most inner 
regions, in cyan). Nonetheless, such a region exists also for a value of the light Higgs mass 
excluded by LEP for the SM, m/j^ = 100 GeV, but only for big values of the mixing angle. 
No new regions (with respect to the SM) in which the model can survive up to the Plank 
scale open for nih^ > 160 GeV, as the allowed space deforms towards smaller values of m/i^. 

D. Heavy neutrino mass influence 

As stated in section IIIIl the RH neutrinos play for the extra scalar singlet the role of the 
top quark for the SM Higgs. This is particularly true for the vacuum stability condition, as 
the fermions in general provide the negative term that can drive the scalar couplings towards 
negative values. Figure [6] shows how the allowed regions in the nih^-mh^ plane change for 
a RH Majorana neutrino Yukawa coupling y^'^ = 0.2 (that for x = 3.5 TeV correspond to 
rriuf^ = 1 TeV), not negligible anymore. For y*^ = 0.4, the changes are even more drastic, 

^ The chosen values maximise the probabihty for the decays hi — >■ bb, hi 77, hi — >■ W'^W' and hi — ^ ZZ, 
respectively. 
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FIG. 4: Allowed values (that are those between the same colour lines) for nih^ as a function 
of the scale Q in the B — L model by eqs. (27^ and ^ME), for several values of ruh^ and (^a[) 
a = 0, ( f^6p a = 0.1, ( [Tcp a = tt/8 and (Jd^ ) a = tt/A. Also, x = 3.5 TeV and rriy^ = 200 
GeV. Only the allowed values by our convention rrih^ > rnt^ are shown. 



shrinking the allowed region even further. 

The effect of having non negligible y'^^ couplings is evident if we compare figure [6] to 
figure [31 Notice that also the VEV x is changed (from 7.5 TeV to 3.5 TeV), but this is only 
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FIG. 5: Allowed values in the vs. a space in the B — L model by eqs. (21) and ^21 
for (E^j = 100 GeV, (E^j nih, = 120 GeV, ^ ruh, = 160 GeV and (E^j rrih, = 180 
GeV. Golours refer to different values of Qj GeV: blue (^lO'^j, red {\^'^ ), green (^10^°^), purple 
and cyan (VS^'^). The plots already encode our convention vtih^ > nih^ and the shaded 
red region refers to the values of a forbidden by LEP. Here: x = 3.5 TeV, rriy^^ = 200 GeV. 



responsible for the smaller upper bounds of rrih^ in figures [6a] and [6bl For small values of 
a it is evident our analogy between the top quark and the RH neutrinos, as now m^^ has 
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FIG. 6: Allowed values in the nih^ vs. rrih^ space by eqs. (21) and [2^] . for [6o\] a = and 
(EBi) a = 0.1, ( fgcj) a = 7r/4 and (6d\) a = n/3, for rriy^ = 1 TeV and x = 3.5 TeV. Colours 
refer to different values of Q / GeV: blue (lOr'), red (10^), green (lOr"^), purple (10^^) and 
cyan (10^'^ ). The shaded black region is forbidden by our convention > rnh^, while the 
shaded red region refers to the values of the scalar masses forbidden by LEP. 



a sensible lower bound too. The analogy holds also for bigger values of the angle, as the 
allowed region of masses is shrunk from below as we increase the RH Majorana neutrino 
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Yukawa coupling, while the upper bound stays unaffected. The effect is even more evident 
for big values of the scale Q, with the Plank scale precluded now for whatever Higgs boson 
masses at a = tt/A and tightly constraining the allowed ones at a = vr/S. 

Moving to the rrih^-a scan at fixed rrih-^ values, figure [7] shows the effect of the heavy 
neutrinos in this case, to be compared to figure O It is evident that this model can survive 
until very large scales Q with massive heavy neutrinos (for which, y'^^ > 0.2) only for the 
light Higgs boson masses allowed in the case of the SM, that is, rrih^ ~ 160 GeV. The mixing 
angle must also be small, a < tt/S, providing a tight constraint on m/ij. For smaller hi 
masses, the effect of a large y^^ is to preclude scales Q > 10^ GeV almost completely, with 
for example just a tiny strip for rrih^ = 120 GeV for which there exists a combination of 171^2 
and a such that the model is consistent up to Q = 10^'^ GeV. Finally, figure [7d]is not visibly 
different from figure [5d] just because we are showing only the m/ij > m/^ region, the shrunk 
region being below. 

E. VEV effect 

The last effect to evaluate comes from changing the values for the B — L breaking VEV 
X. Figure [Ha] shows the allowed regions in the nih^ vs. a plane for fixed nih^ = 160 GeV 
and y*^ = 0.2 (that is, a particular case that shows all the interesting effects at once). As 
expected, since A2 is a function of rrih^/x (see for instance eq. ([H])), at a = the bounds on 
m/i2 simply scale hnearly with the VEV. Regarding the upper bound, increasing the VEV x 
naively increases the allowed region of the heavy Higgs masses, but it is remarkable that the 
effects are present only for small angles, a < 0.1 radians, being the bigger angles unaffected. 
Concerning the lower bound, or the vacuum stability of the model, at fixed y^^ , increasing 
the VEV X requires to increase rrih^ to keep A2 constant at the EW scale. This explains why, 
with non negligible y^ , the allowed heavy Higgs masses are shrinking from below when we 
increase the VEV x, as one can see in figure IHal and comparing figure l8b] with figure [6c| both 
for a = it/A and y^^ = 0.2, but for x = 3.5 and x = 7.5 TeV, respectively. 

In general, for the model to survive up to very large scales Q ~ Mpianck, it is preferred 
the heavy neutrinos to be light with respect to the VEV x, in such a way that their Yukawa 
couplings are negligible in the RGE evolution of the scalar sector. 
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V. SUMMARY AND CONCLUSIONS 



We have investigated the triviahty and vacuum stabihty conditions of the minimal (or 
pure) B — L model with a particular view to define the phenomenologically viable regions of 
the parameter space of the scalar sector, by computing all relevant RGEs (gauge, scalar and 
fermionic) at the one-loop level in presence of all available experimental constraints. The 
RGE dependence on the Higgs masses and couplings (including mixings) has been studied 
in detail for selected heavy neutrino masses and couplings as well as discrete choices of the 
singlet Higgs field VEV. 

Altogether, we have found that there exist configurations of the model for which its 
validity is guaranteed up to energy scales well beyond those reachable at the LHC while at 
the same time enabling the CERN hadron collider to probe its scalar sector in Higgs mass 



22 



and coupling regions completely different from those accessible to the SM. Furthermore, we 
have shown that investigations of the Higgs sector of this extended scenario may also lead to 
constraints on other areas, such as the (heavy) neutrino and Z' sectors (the latter indirectly, 
through the VEV of the singlet Higgs state directly intervening in the scalar RGBs). 

Combining the results of this paper on triviality and vacuum stability with those on 
unitarity of Ref. 18|, we are now in a position to investigate the production and decay 
phenomenology of both Higgs states of the minimal B — L model at present and future 



accelerators 



3l|. 



Appendix A: RGEs 

In this appendix we present the complete set of one-loop RGEs for the minimal U{1)b-l 
extension of the SM. For some parameters, the equations will be equal to those of the SM, 
as no extra contribution arises at one-loop level. 

1. Gauge RGEs 



The RGEs for the SU{?>)c and SU{2)l gauge couplings gs and g are |32| 
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di 
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where rig = 3 is the number of generations. 

Following standard techniques, we obtain for the Abelian couplings 
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[A'^'^91] , 

[A^^g'^ + 2A^^gfg + A^^g[g'] , 

[A^'^gig' + 2gf) + 2A^^ g[{^ + gl) + A^^g[^g] 



with 



(Al) 
(A2) 



(A3) 
(A4) 
(A5) 

(A6) 



where the first sum is over the left-handed two-component fermions and the second one is 
over the complex scalars. For the model we are discussing {Y is the SM weak hypercharge. 



23 



X = B — L is the B — L number), the coefficients of eq. (1A6P are, respectively: 

A^^ = 41/6, A^^ = 12, A^^ = 16/3. (A7) 



2. Fermion RGEs 



From straightforward calculations we obtain: 



d ^ yt f9^2 o 2 2 17 2 17^ 2 ,3 5, 



= _ 8g' - -g' - —gf - —g' - -g,' - -gg[ . (A8) 

dr 167r2 V2^* 4^ 12^^ 12^ 3^^ 3^^V 

For the right-handed neutrinos, it is not restrictive to consider the basis in which the Ma- 
jorana matrix of couplings is real, diagonal and positive: = diag (yf^, y^, y^). Then we 



get |33|,l34| 



7. 

" „M Hi 



dr' 167r2 



(4(yf )2 + 2Tr [{y'^f] - Qg'^) , = 1 . . . 3) . (A9) 



3. Scalar RGEs 



A very straightforward way to find the one-loop RGEs for the parameters of the scalar 
potential is to compute the one-loop effective potential and to impose its independence from 
the renormalisation scale. To one-loop level, the scalar potential V reads: 

V = + AV(^) , (AlO) 

where V^'^^ is the tree-level potential and AV^^^ indicates the one- loop correction to it. To 
compute the latter it is useful to re-write the tree-level potential 

V^''\H, x) = m^H^H + /i' I X I' +Xi{H^Hf + A2 | x I' +X3H^H \ x T (All) 
in terms of the real scalar fields: 

X= ^(05 + ^06). (A12) 




^ Notice the we get a difference of a factor 3 in the third term in the RHS of the last expression in eq. (14) 



contained in Ref. [34|. The authors of Ref. |34| acknowledged the difference and will correct their paper. 
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The only combinations of fields that are involved are (jP' = 0i + 02 + 03 + 04 = 05 + 06' 



so that eq. flAlip becomes: 



(A13) 



The one-loop correction to the tree-level potential flAlSP is, in the Landau gauge, 



• Mf(0^r/^) 
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where c, are constants that depend on the renormalisation scheme (for example, in the MS 
scheme, it is q = 3/2 for scalars and fermions, q = 5/6 for vectors). Expanding eq. (lA14p 
and keeping the contributions of the scalar fields (Higgs and Goldstone bosons), of the 
top-quark, of the gauge bosons and of the RH neutrinos only, we obtain 
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where the field-dependent squared masses are, in a self-explanatory notation: 
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As usual, we define the beta functions Pi (z = 1 . . . 3) for the quartic couplings, the gamma 
functions '~frn,fi for the scalar masses and the scalar anomalous dimensions 7<^_^ as follows 
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(t = lnQ): 
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dii^ 

~dt 

dt 

drf 

~dt 

Now we can extract the RGEs for the parameters of the scalar potential just by requiring 
that the first derivative of the effective potential with respect to the scale t vanishes: 
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dt dt^ ' 



(A26) 



keeping only the one-loop terms. Reorganising it in a more convenient way, we see that 
eq. ( ]A26p implies the following equations: 
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Imposing that each term between squared brackets vanishes, we can obtain the RGEs 
for the parameters of the scalar potential after inserting the explicit expression of the scalar 



anomalous dimensions 7,^ and 7^. The latter are easily computed and read 

1 
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Inserting eqs. ( ]A27p and (1A28P into the RGEs, we finally obtain the RGEs for the five 



parameters in the scalar potential: 

1 dm"^ 1 ^9 ^M^x Qo 3, \ / . 
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